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Abstract 

The leading term in the semileptonic width of heavy flavor hadrons depends on 
the flfth power of the heavy quark mass. We present an analysis where this power 
can be self-consistently treated as a free parameter n and the width can be studied 
in the limit n — oo. The resulting expansion elucidates why the small velocity 
(SV) treatment is relevant for the inclusive semileptonic b ^ c transition. The 
extended SV limit (ESV limit) is introduced. The leading terms in the perturbative 
as expansion enhanced by powers of n are automatically resummed by using the 
low-scale Euclidean mass. The large-n treatment explains why the scales of order 
mh/n are appropriate. On the other hand, the scale cannot be too small since the 
factorially divergent perturbative corrections associated with running of show up. 
Both requirements are met if we use the short-distance mass normalized at a scale 
around rrih/n ~ 1 GeV. A convenient deflnition of such low-scale OPE-compatible 
masses is briefly discussed. 
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1 Introduction 



A true measure of our understanding of heavy flavor decays is provided by our ability 
to extract accurate values for fundamental quantities, like the CKM parameters, 
from the data. The inclusive semileptonic widths of B mesons depend directly on 
\Vcb\ and \ Vub\, with the added bonus that they are amenable to the OPE treatment 

1^, 0, and the non-perturbative corrections to rsi(-B) are of order I /ml |^, |, 0, 
i.e. small. Moreover, they can be expressed, in a model independent way, in terms of 
fundamental parameters of the heavy quark theory - fi^ and fiQ - which are known 
with relatively small uncertainties. 

Therefore, the main problem in the program of a precise determination of, say, 
\Vcb\ from rsi(i?) is the theoretical understanding of the perturbative QCD corrections 
and the heavy quark masses. The theoretical expression for rsi(i?) depends on a 
high power of the quark masses rrif, (and m^): Ts\ oc with n = 5. The quark 
masses are not observables, and the uncertainties in their values are magnified by 
the power n. For this reason it is often thought that the perturbative corrections in 
the inclusive widths may go beyond theoretical control. 

The quark masses in the field theory are not constant but rather depend on the 
scale and, in this respect, are similar to other couplings like a.^ which define the 
theory. Of course, mg and as enter the expansion differently. Moreover, in contrast 
to as in QCD the heavy quark mass ttiq has a finite infrared limit to any order in 
the perturbation theory, mq'^^, which is routinely used in the calculations. 

If rsi(_B) is expressed in terms of the pole mass of the b quark treated as a given 
number, the expression for rsi(i?) contains a factorially divergent series in powers 
of as 



due to the l/mg infrared (IR) renormalon |^, ^, |TD[. This series gives rise to an 
unavoidable uncertainty which is linear rather than quadratic in l/rrib. The related 
perturbative corrections are significant already in low orders. They are directly asso- 
ciated with running of as and their growth reflects the fact that the strong coupling 
evolves into the nonperturbative domain at a low enough scale. A strategy allowing 
one to circumvent this potentially large uncertainty was indicated in Refs. [H, ^ 
- instead of the pole quark mass one should pass to a Euclidean mass, according 
to Wilson's OPE. This mass is peeled off from the IR part; the IR domain is also 
absent from the width up to effects l/rriQ, and the (IR-related) factorial divergence 
disappears. Following this observation, it became routine to express rsi(-B) in terms 
of the Euclidean quark masses in the MS scheme, say, nibirnb). 

Killing IR renormalons on this route, however, one does not necessary have a 
fast convergent perturbative series: in general, it contains corrections of the type 
{nasY which are not related to running of as - they are present even for the van- 
ishing (3 function. We will see that such large corrections inevitably appear if one 




1 



Figure 1: The decay distribution over invariant mass of the lepton pair and over 
energy release E,. = mi, — rric — {(f)^^"^ in b ^ ciu at rric/mh = 0.3. 

works with the MS masses hke Tni){mb) and fnc{mc). Due to the high value of the 
power n they constitutes quite an obvious menace to the precision of the theoretical 
predictions. The standard alternative procedure of treating inclusive heavy quark 
decays, which gradually evolved in the quest for higher accuracy, is thus plagued by 
its own problems. 

The aim of the present work is to turn vices into virtues, by treating Tl clS cL 
free parameter and developing a. 1/n expansion. In this respect our approach is 
conceptually similar to the l/N^ expansion or the expansion in the dimensions of 
the space-time, etc. which are quite common in various applications of field theory 
(for a discussion of the virtues of expansion in 'artificial' parameters, see [llj). 



The main advantage is the emergence of a qualitative picture which guides the 
theoretical estimates in the absence of much more sophisticated explicit higher-order 
calculations. 

Certain unnaturalness of the MS masses normalized, say, at /i = m^, for inclusive 
widths is rather obvious before a dedicated analysis, in particular in b ^ ciu inclu- 
sive decays. The maximal energy fed into the final hadronic system and available for 
exciting the final hadronic states, which determines the "hardness" of the process, 
is limited by rrih — rric — 3.5 GeV. Moreover, in a typical decay event, leptons carry 
away a significant energy Eei, > y/q^, since the lepton phase space emphasizes the 
larger g^. This is illustrated in Fig. 1 where the distribution over invariant mass of 
the lepton pair and energy release Ej- = rUb — rric — y/(f is shown. The situation 
is less obvious a priori for b ^ u, but again the typical energy of the final state 
hadrons is manifestly smaller than m^. We use the parameter n to quantify these 
intuitive observations. 

Since the power n = 5 is of a purely kinematic origin, technically it is quite easy 
to make n a free parameter. To this end it is sufficient, for instance, to modify the 
lepton current by introducing fictitious additional leptons emitted by the W boson, 
along with the standard iu pair. At the very end the number of these fictitious 
leptons is to be put to zero. 
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The emergence of the large perturbative terms containing powers of n in a general 
calculation is rather obvious and can be illustrated in the following way. Consider, 
for example, the 6 — > w decay rate, and use the pole mass The series will have 

a factorial divergence in high orders due to the 1/m IR renormalon, but we are not 
concerned about this fact now since we reside in purely perturbative domain and we 
do not study high orders /c ~ l/a^ of the perturbation theory. We can consider, for 
example, the case of vanishing /? function when nothing prevents one from defining 
the pole mass with the arbitrary precision. 

Let us assume then that using the pole mass, the perturbative expansion 

of the width 
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has all coefficients completely n-independent (we will show below that it is the 
case in the leading-n approximation). Then, being expressed in terms of a different 
mass TTib, 

(2) 
(3) 
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Just using a different mass generates n-enhanced terms ~ {na^^ . In the case of 
the MS mass one has, basically, c — 4/3. In order to have a good control over the 
perturbative corrections in the actual width, one needs to resum these terms (at 
least, partially). This can be readily done. 

These n-enhanced perturbative corrections are only one, purely perturbative, 
aspect of a general large-n picture. The total width is determined by integrating n- 
independent hadronic structure functions with the n-dependent kinematic factors, 
the latter being saturated at large \fq^ close to the energy release my, — rric- As 
a result, the energy scale defining the effective width of integration, is essentially 
smaller than the energy release; a new, lower momentum scale automatically emerges 
in the problem at large n. 

The question of the characteristic scale in inclusive decays was discussed in the 
literature. The discussion was focused, however, almost exclusively on the problem 
of choosing the normalization point for the running coupling a^. We emphasize 
that for the inclusive decays it is a secondary question; the principal one is the 
normalization point for quark masses which run as well. Although their relative 
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variation is smaller, this effect is enhanced, in particular, by the fifth power occurring 
in the width. 

With n = 5 treated as a free parameter, we arrive at the following results con- 
cerning semileptonic B decays driven by 6 ^ luq, q = c or u: 



1. The leading n-enhanced perturbative corrections to the total width are readily 
resummed by using the low-scale quark masses. No significant uncertainty in 
the perturbative corrections is left in the widths. 

2. The surprising proximity of the b —>■ ciu transition to the small velocity (SV) 



limit [|T^] which seems to hold in spite of the fact that m^/m^ -C 1, becomes 
understood. The typical kinematics of the transition are governed by the 
parameter {'mh — mc)/nmc rather than by {nib — irtc) / mc- This "extended SV" 
parameter shows why the SV expansion is relevant for actual b ^ c decays. 

3. The 'large-n' remarks are relevant even if as does not run. In actual QCD, 
the better control of the perturbative effects meets a conflict of interests: the 
normalization point for masses cannot be taken either very low or too high, 
thus neither MS nor pole masses are suitable. We discuss a proper way to 
define a short-distance heavy quark mass mgi^fi) with /i <^ rriQ which is a 
must for nonrelativistic expansion in QCD. A similar definition of the effective 
higher- dimension operators is also given. 

4. The limits mq oo and n —>■ oo cannot freely be interchanged. Quite different 
physical situations arise in the two CclSGS, clS indicated later. The perturbative 
regime takes place if {nib — rnc)/n ^ Aqcd, and the process is not truly short- 
distance if the opposite holds. Nevertheless, at nic ^ Aqcd even in the latter 
limit the width is determined perturbatively up to effects ~ AQQj^/ml. 

We hasten to add that analyzing only terms leading in n may not be fully ade- 
quate for sufficiently accurate predictions in the real world where n = 5, which turns 
out to be not large enough. Still, it is advantageous to start from the large n limit 
and subsequently include (some of the) subleading terms. 



2 The theoretical framework: large-n expansion 
in inclusive semileptonic decays 

The new tool we bring to bear here is the following: the actual value n = 5 of the 
power n in the width arises largely through the integration over the phase space 
for the lepton pair. One can draw up simple physical scenarios where n appears 
as a free parameter and the limit n — *• oo can be analyzed. The simplest example 
is provided by / (massless) scalar "leptons" 0, emitted in the weak vertex of the 
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lepton-hadron interaction, 



G 

where Gi generalizes the Fermi couphng constant to / 7^ 0. At / = the couphng 
Go = Gp- Then by dimensional counting one concludes that T{Q) ~ \Gi\'^ ■ ttiq, 
where n = 2/ + 5. It should be kept in mind that the details of the physical scenario 
are not essential for our conclusions. 

We will discuss b ^ q transitions with an arbitrary mass ratio, mg/mf,, thus 



incorporating both h ^ u and b ^ c decays. Following Refs. |0, |T3] we introduce 
a hadronic tensor h^^i^qQ, g^), its absorptive part Wfj^y^qo, q^) = 4 disc hfj_ty{qo, q^) and 
its decomposition into five covariants with structure functions Wi{qo, q'^), i = 1, 5. 
Variables go and are the energy and the effective mass of lepton pair. Only wi and 
W2 contribute when lepton masses are neglected and the semileptonic decay width 
is then given by 



where 



^^2nJo ^%y%-(f\(lWi{qo,q) + -{qQ-q)w2{qo,q)^, (6) 

9max = {Mb - Md) , gomax = ^J^^ • 

Note that the upper limits of integration over g^ and go are determined by vanishing 
of the structure functions for the invariant mass of the hadronic system less than 
M|). The lower limit of integration over go and g^ is due to the constraints on 
the momentum of the lepton pair. It has nothing to do with the properties of the 
structure functions. In particular, the structure functions exist also in the scattering 
channel where the constraint on the lepton momentum would be different. In what 
follows we will see that the distinction between the leptonic and hadronic kinematical 
constraints is important. 

By adding I extra scalar "leptons" emitted in the weak vertex we make n = 5 + 2/ 
a free parameter. The quantity 7 then becomes n dependent, 7 — > '^{n), 

1 r^max ^ ^ , /'QOmax I ~ T 9 9 

= — dq' {q'y / dgo ^gg - g^ {g'wi(go, g')+ 
Ztt Jo J\ ^ 



^{ql-q')w2{qo,q') \ ■ (7) 



3 

In the real world / = 0, = 5, but as far as the QCD part is concerned we are free 
to consider any value of I. The extra factor (g^)' appeared due to the phase space 
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of the scalar "leptons" .[] 

A different choice of variables is more convenient for our purposes. Instead of go 
and we will use e and T, 



e = Ms - go - ^Ml + ql - q\ T= ^ Ml + gg - q^ - Mj, . (8) 

The variable T is simply related to the spatial momentum q. More exactly, T has 
the meaning of the minimal kinetic energy of the hadronic system for the given 
value of g. This minimum is achieved when the D meson is produced. The variable 
e has the meaning of the excitation energy, e = {Mx — Md) + {Tx — T) where 



Tx = \J Mjr + qQ — q^ — Mx is the kinetic energy of the excited state with mass 
Mx- The following notations are consistently used below: 

A = Mb - Mfl, g' = go - g' , |g | = 
In terms of the new variables one obtains 

= - / dT (T+Md) + 2MdT / de ( - 2MbT - 2Ae + 2Te + e 
n Jo Jo ^ 

|(A2 - 2MbT - 2Ae + 2Te + e^)wi + ^(T^ + 2MdT)w2^ (9) 

where 

A2 / 

T„.ax = TTirr ' e^ax = A - T - ^T^ + 2M^T ; 
zivib 

it is implied that the structure functions Wi^2 depend on e and T. 

2.1 Cancelation of the infrared contribution 

Before submerging in the large n limit we will address a question which naturally 
comes to one's mind immediately upon inspection of Eq. (|^) or Eq. (|^. Indeed, 
these expressions give the total decay probabilities in terms of an integral over the 
physical spectral densities over the physical phase space. Both factors depend on 
the meson masses, and know nothing about the quark mass. This is especially 
clear in the case of the phase space, which seems to carry the main dependence for 
large n. And yet, in the total probability the dependence on the meson masses must 
disappear, and the heavy quark mass must emerge as a relevant parameter. In other 
words, the large-distance contributions responsible for making the meson mass out 
of the quark one must cancel each other. 

^ Strictly speaking, the extra "leptons" produce a change in the leptonic tensor, in particu- 
lar, making it non-transversal. This variation leads, in turn, to a different relative weight of the 
structure function W2 , and the emergence of the structure functions Wi and in the total proba- 
bility. These complications are irrelevant, and we omit them. We also omit an overall I dependent 
numerical factor, which merely redefines G;. 
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It is instructive to trace how this cancelation occurs. We will study the issue 
for arbitrary value of I, assuming for simplicity that the final quark mass vanishes 
(i.e. h ^ u transition). This assumption is not crucial; one can consider arbitrary 
ratio mc/mi), with similar conclusions. Since the main purpose of this section is 
methodical we will limit our analysis to terms linear in A. The basic theoretical 
tool allowing one to trace how in the decay probability is substituted by m^, 
as a result of cancelation of the infrared contributions residing in Mb and in the 
integrals over the spectral densities, is the heavy flavor sum rules presented in great 
detail in Ref. IH, see Eqs. (130) - (135). 

We start from expanding the integrand in Eq. (^ in e, keeping the terms of the 
zeroth and first order in e, 



7(n) 



- / dT T^Mjj (Mb -2r)'x 

TT Jo 



Mb{Mb-2T) I dewi-2{l + l){MB-T) I deewi 







+ 



1 fMB/2 

— dTT^M^B\MB-2Ty-^x 
3tt Jo 

pMb-2T i-Mb~2T 

Mb{Mb-2T) dew2-2l{MB -T) deewa 

Jo Jo 



The integrals over e are given by the sum rules mentioned above, 
■^Jdewi = l, ^ J deewi = Mb - rUb = A , 

and 



1 



dew2 



2mh 



1 



deew2 



2my 



A. 



(10) 
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271 J T 27T J T 

The expressions (p^ ) and (0) imply that the integration over e is saturated at small 

e, of order A. After substituting the sum rules in Eq. ( p!0D and integrating over T 

we arrive at _ 

A 



7(n) 



2(/ + 4)(/ + 3)(/ + 2) 



1 



{21 + 5) 
{21 + 4) 



Mb 
A 



M 



B 



(13) 



2(/ + 4)(/ + 3)(/ + 2)(Z + l) 

The first term here comes from Wi and the second from W2- This expression explicitly 
demonstrates that the meson mass is substituted by the quark one at the level of 
1/m terms under consideration, Mb{1 — A/ Mb) = rrib. 

Thus, the cancelation of the infrared contribution is evident. It is worth em- 
phasizing that the infrared contributions we speak of need not necessarily be of 
nonperturbative nature. The cancelation of infrared l/mj, effects takes place for 
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perturbative contributions as well as long as they correspond to sufficiently small 
momenta, ^rrib/n. 

Let us parenthetically note that Eq. (^) is valid for any /; in particular, we 
can put I equal to zero. One may wonder how one can get in this case a non- 
vanishing correction associated with / dee 1/72, while the weight factor for W2 in the 
original integrand at Z = seemingly has no dependence on e at all. This case 
is actually singular: the domain of T — > Mb/2 contributes to the integral at the 
level Aqcd/ Mb- For such T the interval of integration over e shrinks to zero, but 
nevertheless the integral of W2 stays finite (unity) according to Eq. (|l^); this signals 
a singularity. It is directly seen from Eq. ( [T0| ) where at / = one has logarithmic 
divergence at T — Mb/2. The most simple correct way to deal with at / = is 
inserting a step-function 9 {Mb — 2T — e) in the integrand instead of the upper limits 
of integration. This step-function does provide an e dependence, the term linear in 
e is e6{MB — 2T). At non- vanishing values of / the 5-function above produces no 
effect. As usual, analytic continuation from the nonsingular case / > to / = leads 
to the same result. 

In principle, it is instructive to trace the / dependence of other non-perturbative 
terms, e.g. those proportional to fi^ and fi^ which are known for 6 — > m at arbitrary 
I from explicit calculations of Refs. §|, and have a very simple form. We will 
not dwell on this issue here. 

Our prime concern in this work is the interplay of the infrared effects in the 
perturbative corrections. The main new element appearing in the analysis of the 
perturbative corrections is the fact that the integrals over e are not saturated in the 
domain e ~ A. That is why we need to introduce the normalization point fi which 
will divide the entire range of the e integration in two domains, e > /i and e < /i. 

The domain e < fi provides us with a clear-cut physical definition of A, which 
becomes dependent, 

+ 21,-1+ 3|f| y m.j^^yq') J„"d6»f-"(f,lfl) ' * ' 

where the variable e is defined by Eq. (||) with Md set equal to zero, the second 
variable T = \q \ in the h ^ u transition. Equation (|1^ is a consequence of the sum 



rules derived in Ref. [|14|. There, the integrals over the spectral densities are given in 
the form of a "condensate" expansion. The currents inducing given spectral densities 
(and the spectral densities themselves) acquire certain normalization factors due 
to short-distance renormalization of the weak vertex, which are irrelevant for our 
present purposes. To get rid of these normalization factors we consider the ratio of 
the sum rules. The power corrections to the denominator in the right-hand side of 
Eq. ([1^ shows up only in O (Aqcp/g^^ terms. 

The quark mass, mf,{jj,) = Mb — A(yu), becomes in this way a well-defined and 
experimentally measurable quantity. A very similar definition of A(/i) does exist 
also in the b ^ c transition, through the integral over w\~^^, see Ref. This 
definition generalizes the Voloshin sum rule [l^ to the arbitrary values of q, and 
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includes corrections O^K^qy)/'''^)- See Section 4 for the further discussion of A(yu) 
definition. 

2.2 The large n limit 

Now when we understand how the infrared parts cancel, we are ready to address 
the practical issue of what particular normalization point n is convenient to use in 
the analysis of the total widths. Certainly, the theoretical predictions can be given 
for any value oi Aqcd- Our goal is to choose /i in such a way that the domain 
of momenta above /i does not give the enhanced perturbative corrections. Then all 
perturbative contributions enhanced by large n will be automatically included in 
the definition of m{,(/i). 

To this end we invoke the large n limit, as was explained in the Introduction. 
For / ^ 1 one has 



( - 2MbT - 2Ae + 2Te + e^)' ~ A^' exp 



-21 



(TMu e Te 



\^ A2 A A2 2A2 
and the dominant domains in the integration variables are given by 

All expressions then simplify, and the two integrations decouple 



7(n) 



Stt Ml 



^"°drr^r + 2ro)Mr + ro)e-('+^)^/^°°dee-2(^+^)i^i[e, ^) + 



A^ f°° 3 3 ( tA'^ \ 1 

drr 2 (r + 2ro)^ (r + tq) e^'" / de q-'^^^w^ e, — — (15) 



l2MlJo ■ ' ' '"'^ Jo " ^y'2MB, 

where 

_ 2MbT _ 2MdMb 
The integrals over the excitation energy e 



1 

J,{T- a) = — dee-^u7,(e, T) (16) 



27T Jo 



are a combination of the sum rule considered in Refs. [0, They are related to 
the Borel transform of the forward scattering amplitude of the weak current off the 
B meson. The quantity a is the Borel parameter, and expanding Jj(T; a) in 1/a 
gives us the first, second, and so on sum rules of Ref. |T^. The Borelized version 



has an advantage, however, of providing an upper cut off in e in a natural way. This 
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is important in the analysis of the perturbative corrections; the second argument of 
J defines the normahzation point. 

One readily reads off from Eq. (ITH) that only a fraction ~ 1/n of the total energy 



release is fed into the final hadronic system. The three-momentum carried by the 



final state hadrons is likewise small, namely of order ymcA/n for 6 ^ c and A/n for 
b ^ u, respectively; the latter case is obtained from the former by setting M^, = 0. 
This implies that by evaluating the various quantities at a scale ~ A/n rather than 
A one includes the potentially large higher order corrections. These conclusions can 
be made more transparent by considering two limiting cases. 

(a) If rric stays fixed, the limit n ^ oo leads to the SV regime, as is expressed 
by w ~ {nib — mc)/nmc <^ 1. Neglecting r compared to tq, one gets 



--it)^'^r-— -•®^)- 



7(n) 



3 \iviB/ "'U \z,iviQ 

The behavior of Fsi at n — > oo then is determined by the integrals Ji^2 near zero ar- 
gument. The second term is clearly subleading because the weight function contains 
an extra power of r. It is not difficult to see that the axial current contribution is 
dominant, and one obtains for the reduced width 

0»|_.A"v^(^) (18) 

where ^A{f^) = 1 + 0(Q;s(/i)) is the coefficient function of the unit operator in the 
first sum rule for the axial current at zero recoil (for further details see Ref. fl^). 
To the degree of accuracy pursued here (power corrections are switched off so far) 
one has C,a — Va where the factor rjA incorporates the perturbative corrections to 
the axial current at zero recoil.| Eq. (|18D provides a decent approximation to the 



true width even for n = 5. The vector current contribution is subleading; to take it 
into account at n = 5 we need a refined expansion to be described in Appendix. 

We note that the leading-?7, expansion of the width yields {Mb — M^)" which 
differs from {mh—rric)^ by terms {n Og)^ that we target, but coincides to this accuracy 
with {mb{fi) — mc(/i))" if fi ^ rrih/n, 

mt,{fi) - m,{fi) ^ {Mb - M^) + (m, - ^J^flMlL^M + ... 

Zmbrric 

c^{MB-Mn) fl-^Tr^+ ■■■] ■ (19) 



■^The quantity r]\ = lim^^o Ca'''^*(/^) is ill-defined once power corrections are addressed; at the 
same time ^a(/^) is even then a well defined quantity provided that ^ 3> Aqcd- Yet, to any finite 
order in perturbation theory, one can work with rjA- 
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In the last line we used the perturbative expression for — /^^(/i)- 

(b) Another case of interest is rUg <^ rrib/n, when the final-state quark is ultrarel- 
ativistic. In this case the axial and vector currents contribute equally. The reduced 
width now takes the form 

.(.„,^iMB/:d..^e-...(l|2;^). (.0, 

The quantity Ji{tMb/2) stays finite (unity at the tree level) at r — 0. The second 
structure function again formally yields only subleading in 1/n contributions. Thus 

7n(6^«)|™ ^ M'^^UMb/u), (21) 

where 

un) = J. ^) , er = 1 ■ (22) 

\ n n J 

This width decreases faster with n than for massive quarks due to the higher power 
of |g I ~ m^t/n. This underlies the fact that the numerical factor in Fgi in front of 
is much smaller for b ^ u than in the SV limit, namely, 1/192 vs. 1/15. For the 
same reason the simple expansion described above provides a poor approximation 
for 6 ^ M when evaluated at Z = 0, i.e. n = 5. 

(c) One can consider the third large-n regime when mc/rrif, ~ 1/n. Although it 
appears to be most close to the actual situation, we do not dwell on it here: the 
analysis goes in the very same way, but expressions are more cumbersome since one 
has to use the full relativistic expression for the kinetic energy of the D meson. No 
specific new elements appear in this case. 

(d) Let us discuss now the normalization point for quark masses at large n. 
Eqs. (|18|), ( |22D show the kinematically-generated dependence on the masses. The 
non-trivial hadronic dynamics are encoded in the factors and ^u- Our focus is the 
heavy quark masses. In the SV case (a) the dependence on them is trivial as long as 
the low-scale masses are employed. In the case of the 6 — ^ -u transitions we showed 
in Sect. 2.1 that the factor ^u{MB/n) actually converted into mfc(/i)", and, thus, 
effectively established the normalization point for the latter, /z ^ rrih/n. Using the 
effective running mass mfe(/i ~ mf,/n) does not incorporate, however, the domain of 
the gluon momenta above ~ rrih/n. This contribution has to be explicitly included 
in the perturbative corrections. It is not enhanced by powers of n. Physically, it 
is nothing but a statement that for r{b u) the proper normalization scale for 
masses is nib/n. 

We hasten to emphasize that the statement above is not a question of normal- 
ization of as used in the perturbative calculations. One can use at any scale 
to evaluate mi,{fi) as long as this computation has enough accuracy. The fact that 
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using inappropriate can lead to an apparent instability in mh{fi), is foreign to the 
evaluation of the width. 

Let us summarize the main features which are inferred for inclusive widths by 
analyzing the straightforward large-n expansion introduced via Eq. (]TB|). The in- 
tegral over the lepton phase space carries the main dependence on n. The QCD 
corrections - the real theoretical challenge - are only indirectly sensitive to n: for 
the kinematics determines the energy and momentum scales at which the hadronic 
part has to be evaluated. 

The characteristic momentum scale /i for the inclusive decay width is smaller 
than the naive guess fi ~ rrib. In the large n limit it scales like rrih/n. This momentum 
defines the relevant domain of integration of the structure function in qq, i.e. the 
scale at which the forward transition amplitude appears in the total decay rate. In 
particular, by evaluating the quark masses that determine the phase space at this 
scale ~ rrib/n one eliminates the strongest dependence of the radiative corrections 
on n. 

The expansion in n derived from Eq. (|1^) allows a transparent discussion of the 
underlying physics. Unfortunately it yields, as already stated, a decent numerical 
approximation only for very large values of n; i.e., for n = 5 the non-leading terms are 
still significant. Not all those terms are dominated by kinematical effects, and their 
treatment poses non-trivial problems. A more refined expansion can be developed 
that effectively includes the kinematics-related subleading contributions and leads 
to a good approximation already for n = 5. This treatment, however, is rather 
cumbersome and less transparent. It will be briefly described in the Appendix. The 
purpose is only to demonstrate that the essential features of our expansion leading 
to the qualitative picture we rely on - that the characteristic momentum scale is 
essentially lower than nif, - hold already for the actual case n = 5. The reader 
who is ready to accept this assertion, can skip the technicalities of the refined 1/n 
expansion in the Appendix. 

3 Applications 

We now briefly consider the consequences of the 1/n expansion for a few problems 
of interest. 

3.1 The extended SV limit 

In many respects one observes that the inclusive b ^ c decays seem to lie relatively 
close to the SV limit. Most generally, the various characteristics of the decay depend 
on the ratio m^/m^; although it is rather small, the actual characteristics often do 
not differ much from the case where it approaches 1, when the average velocity of 
the final-state c-quark was small. The proximity of the inclusive 6 — > c decays to 
the SV limit has two aspects. First, it is obvious that the nonperturbative effects 
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work to suppress the effective velocity of the final state hadrons; the most obvious 
changes are kinematical replacement ruf, Mb and rric — > iVfj)^^)*^^)** when passing 
from quarks to actual hadrons. The impact of increasing the mass is much more 
effective for the final-state charm than for the initial-state beauty. This decreases 
the velocities of the final state hadrons significantly at A = rrib — rric — 3.5 GeV. Yet 
this simple effect would not be numerically large enough were it not considerably 
enhanced by the properties of the lepton phase space, which can be easily seen 
comparing it, for example, with the one for semileptonic decays at fixed = 0. 

On the other hand, without any nonperturbative effects, theoretical expressions 
at the purely parton level are known to favor the SV kinematics even for actual 
quark masses. The simplest illustration is provided by the tree-level phase space. 



1 4 2 fi ,, 

vnt m„ , m: mi m" 



zo{m,, = m,^ 1 - 8^ - 12^ bg ^ + 8^ - ^ . (23) 
Y ml ml ml m^ ml j 

It is most instructive to analyze the sensitivity of this expression to and A = 
mh — mc rather than mt and m-c, as expressed through 

Adzo{mb,A) mbdzo{mb,A) 
= Vo OA ' = ^ dmb ' ^''^ 

where 

Ka + f^b = 5 ■ 

In the light quark limit - m^/m^ — >• - one has 2;o|iight = rnt — O (j^^^ and, 
thus 

/^Aliight = , Kfoliight = 5 . (25) 
In the SV limit, on the other hand, one finds 

64 

zolsv ^ —{rrib-mc)^ . (26) 
5 

Therefore, the SV limit is characterized by 

/^aIsv = 5 , Kfelsv = (27) 

For actual quark mass values - mc/mb ~ 0.28 - one finds 

ka ~ 3 , Kb ^2 (28) 

i.e., even for m^/ml ~ 0.08 ^ 1 one is still closer to the SV than the light-quark 
limit. The "half-way" point - = = 2.5 - lies at m^/m^ ~ 0.05! This is 
a consequence of the large parameter n = 5. A similar pattern persists also for 
the low-order perturbative corrections [l^, |T^, and on the level of the power 
corrections [0, |§, 0]. 



The relevance of the SV approximation in the actual inclusive decays thus finds 
a rational explanation in the 1/n expansion. 



13 



3.2 Resummation of large perturbative corrections 

In computing the perturbative corrections for large n one can potentially encounter 
large n-dependent corrections of the type (nas/vr)^, which makes them sizable for 
the actual case n = 5. The related practical concern of certain numerical ambiguity 



of the one-loop calculations of the widths was raised in Ref. pO|. Fortunately, the 
leading subseries of these terms can readily be summed up. The summation essen- 
tially amounts to using the quark masses normalized at a low scale, yU ~ (rrif,, A)/n. 
As long as one does not go beyond a relatively low order in as one can ignore renor- 
malon divergences and simply use the "/c-th order pole mass" . Quite obviously, the 
terms ~ (nas/vr)'^ appear only if one uses deep Euclidean masses. If the result is ex- 
pressed in terms of low-scale Euclidean masses (i.e. those normalized at yU ~ A/n or 
/i = several units xAqcd) the terms ~ {nag/'n)^ enter with coefficients proportional 
to /i/A and, therefore, contribute only on the sub leading level. 

The above assertion is most easily inferred by applying Eq. (|15D to Wi computed 
through order k (with k < n); they are self- manifest in Eqs. (p^ ) and (p2D. Upon 
using the thresholds in the perturbative Wi to define rrify and rric, the moments of Wi 
are smooth and independent of n in the scaling limit m, A ~ n, which ensures the 
absence of the leading power of n. The perturbative thresholds, on the other hand, 
are given just by the pole masses as they appear in the perturbation theory to the 
considered order. 

It is worth clarifying why these terms emerge if one uses different masses. The 
reason is that the moments of the structure functions (or more general weighted 
averages in the 'refined' expansion) intrinsically contain the reference point for the 
energy go and g^ax- > ?max the integral over go merely vanishes whereas for 

g^ < g^a^x it is unity plus corrections. A similar qualification applies to evaluation of 
the integral over the energy go itself. These properties which single out the proper 
'on-sheir masses, were tacitly assumed in carrying out the expansion in 1/n. 

On the formal side, if one uses a mass other than the perturbative pole mass, 
the perturbative Wi contain singular terms of the type 

S(e - 6m) = 6(e) + V (-6m)''--P- . (29) 

where 6m ~ ■ A is the residual shift in the energy release. (A similar shift in 
the argument of the step-functions emerges for the perturbative continuum contri- 
butions.) Each derivative of the (5-function generates, for example, a power of n in 
the integral over e in the representation ([T5|). For b —>■ u widths considered in detail 

in Sect. 2, using mb{mb) would yield —6m ~ leading to the series (^ff^) /k\ 

which sums into e"^"^/^'^ and converts fni,{mb)^ into m^{^) with -C mb- 

A more detailed analysis can be based on the refined 1/n expansion described 
in the Appendix, Eqs. (^)-(|i5|); it will be presented elsewhere. For our purposes 
here it is enough to note that using masses normalized at a scale ~ A/n does not 
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generate terms ~ (ncts/vr)'^. 

A qualifying comment is in order. Including higher order terms in the expansions 
in 1/n and in a.,, one cannot find a universal scale most appropriate for aZ/perturba- 
tive corrections simultaneously, including the normalization of the strong coupling. 
Even upon resumming all n-dependent effects one would end up with the ordinary 
corrections to the structure functions, e.g. those describing the short- distance renor- 
malization of the vector and axial currents. The relevant scale for such effects is 
obviously ~ nib', we are not concerned about them here since they are not enhanced 
by the factor n and, thus, not large. 



3.3 Power corrections and duality in inclusive widths 

Using n as an expansion parameter one can estimate the importance of the higher 
order nonperturbative corrections by summing up the leading terms in n. The simple 
estimates show, however, that these effects do not exceed a percent level, and thus 
are not of practical interest .0 It is still instructive to mention a qualitative feature 
drawn from this analysis. 

The characteristic momentum scale decreases for large n; at A/n ~ Aqcd one is 
not in the short-distance regime anymore. Correspondingly, the overall theoretical 
accuracy of the calculations naturally seems to deteriorate with increasing n (and/or 
increasing rric up to rrib). However, using the sum rules derived in it is not dif- 
ficult to show that even in the limit n — cx) the width is defined by short-distance 
dynamics unambiguously up to ^/ttlq terms. Moreover, in this limit a straightfor- 
ward resummation of the leading nonperturbative effects yields the width in terms 
of the zero-recoil formfactor |-Fd*P and the physical masses Mb and Mj^^,), 

without additional uncertainties. Thus, in the worst limit for inclusive calculations, 
the theoretical computation of the width merely reduces to the calculation of the 
exclusive B D* transition, the best place for application of the heavy quark 
symmetry 
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The large-n arguments may seem to suggest that the width can depend on the 
meson masses rather than on the quark ones, in contradiction to the OPE. In fact, 
since in the 6 — c transitions one encounters the SV limit when n — oo, the leading 
term is given by {Mb — Mjj)" which differs from (m^ — mc)" only by O (^Aq^j^/mQ^ , 
in accordance with the OPE. The problem can emerge in the next-to-leading order 
in 1/n. These corrections can be readily written up in terms of the transition 
formfactors at small velocity transfer, and take the form of the sums of the transition 



probabilities which one encounters in the small velocity sum rules W^- The latter 



ensure that the terms ~ Aqcd/'^q are absent. By the same token, using the third 

^It is important that here we consider large n for a given order in the perturbative expansion. 

^The b —> cTi/ decays represent a special case: due to the sizable r mass the energy release 
is smaller than in 6 —> c/ii^ and b cev and a SV scenario is realized in a rather manifest way. 
The n-enhanced nonperturbative corrections can then reach the ten percent level, but are readily 
resummed using the exact meson masses in the phase space factors. 
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sum rule, one observes that the leading-n term, nfi'^/rriQ, is absent from the width 
although it is obviously present in {Mb — Md)^ (we assume that the quark masses 
are fixed as the input parameters). Similar properties hold to all orders in 1/n, 
however they are not very obvious in this expansion. 

If the infrared part of the quark masses cancels in the widths, why is it not the 
short- distance masses like MS that emerge? The answer is clear: the excited final 
states in the width enter with the weight ~ exp(— ne/A), and the Coulomb part of 
the mass cancels in the width only up to momenta ~ A/n. 

The above consideration only interprets the known OPE results. The general 
question at which scale duality is violated is more instructive. We immediately see 
that, generically, the characteristic momenta are given by the scale mg/n rather than 
by TTiQ. For example, considering the calculation of the relevant forward transition 
amplitude for b —>■ u in coordinate space |§, ^ , we have 

T ~ (30) 

Its Fourier transform ~ p"logj9^ ( + terms analytic in j9^) is saturated at |x| ~ Tq ~ 
p/n, corresponding to a time separation to ~ n/mq. In the case of the massive final 
state quarks one has to ~ n/A. 

In the semileptonic b c decays one may naively conclude that the scale is 
too low. However, here the heavy quark symmetry ensures vanishing of all leading 
corrections even at A — 0, i.e. in the SV limit; this applies as well to the ESV limit. 
The nontrivial corrections appear only at the 1/A^ or n/{mA) level and they are 
still small. 

The situation is, in principle, different in the nonleptonic decays, where the color 
fiow can be twisted (instead of the color fiow from b to c, it can fiow from b to d). 
In the nonleptonic b ^ cud decay, with A ~ 3.5 GeV, one is only marginally in the 
asymptotic freedom domain and a priori significant nonperturbative corrections as 
well as noticeable violations of duality could be expected. In this respect it is more 
surprising that so far no prominent effects have been identified on the theoretical 
side [^.Q We note that a particular model for the violation of duality in the 
inclusive decays suggested in Ref. |^ explicitly obeys the large-n scaling above. 



The importance of the duality- violating terms is governed by the ratio A/ra, 

T^duaLvioL _ ^_A.p(l-^ln^) ^ ^^^^ 



tree 



with p being a hadronic scale parameter. The corrections blow up when A/n ~ /ihad- 
We hasten to remind, however, that this consideration is only qualitative, since we 
do not distinguish between, say n = 5 and n — 2 = 3. The arguments based on the 
refined expansion in the ESV case suggest that the momentum scale is governed by 
an effectively larger energy than literally given by A/n with n = 5. 

priori the effect of the chromomagnetic interaction in order I/toq could have been significant, 
but it gets suppressed due to the specific chiral and color structure of weak decay vertices [Q, 
properties that are external to QCD. 
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4 Low- scale heavy quark mass; renormalized ef- 
fective operators 

Addressing nonperturbative power corrections in b decays via Wilson's OPE one 
must be able to define an effective nonrelativistic theory of a heavy quark nor- 
malized at a scale fi <^ mg; the heavy quark masses in this theory (which enter, 
for example, equations of motion of the Q fields) are low-scale (Euclidean) masses 
mQ{fi). Moreover, even purely perturbative calculations benefit from using them, 
as was shown above: it allows one to solve two problems simultaneously - resum 
the leading n-enhanced terms and avoid large but irrelevant IR-renormalon-related 
corrections. For the perturbative calculations per se, therefore, the question arises 
of what is the proper and convenient definition of a low-scale Euclidean mass? Of 
course, one can always define it through an off-shell quark propagator. Generically 
this will lead to a mass parameter which is not gauge invariant; this feature, though 
not being a stumbling block, is often viewed as an inconvenience. 

Usually, one works with the mass defined in the MS scheme, which is compu- 
tationally convenient. However, due to the well-known unphysical features of this 
scheme it is totally meaningless at the low normalization point fi = several units 
xAqcd- For example, if we formally write the MS mass at this point 

mf (/.) ^ mfiruQ) (l + ^ In ^) (32) 

we get large logarithms In rrib / fi which do not reflect any physics and are not there 
under any reasonable definition of the heavy quark mass. 

Thus, one must find another definition. We suggest the most physical one mo- 
tivated by the spirit of the Wilson's OPE, i.e. the mass which would be seen as a 
'bare' quark mass in the completely defined effective low-energy theory. 



The basic idea is to follow the way suggested in Ref. |T^, and is based on the 
consideration of the heavy flavor transitions in the SV kinematics, when only the 
leading term in the recoil velocity, ~ v'^, is kept. Note that in Subsection 2.1 we 
have used a similar purpose a different kinematical limit when the final particle is 
ultrarelativistic. 

Let us consider the heavy-quark transition Q ^ Q' induced by some current 
Jqq' = Q'^Q in the limit when rnq qi oo and v = q/rriQi ^ 1 is fixed. The 
Lorentz structure F of the current can be arbitrary with the only condition that 
Jqq/ has a nonvanishing tree- level nonrelativistic limit at v = 0. The simplest 
choice is also to consider Q' = Q (which, for brevity of notation, is assumed in what 
follows). Then one has 



A(^) = lim iMHo-rriQifi)] = lim lim A iL^il^f^lflZ!^ , (33) 
where w{e, \q\) is the structure function of the hadron (the probability to excite the 
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state with the mass + e for a given spatial momentum transfer g ); Mhq is the 
mass of the initial hadron. 

The OPE and the factorization ensures that 

a) For the given initial hadron the above A(/i) does not depend on the chosen 
weak current. 

b) The value mQ(/i) = Mhq — A(/i) does not depend on the choice of the hadron 
Hq. This holds if is taken above the onset of duality. 

To define perturbatively the running heavy quark mass one merely considers 
Eq. (|33D in the perturbation theory, to a given order: 

mQ(/i) = \Mhq] ^ - [Ki^)] ^ . (34) 



pert 



pert 



In the perturbation theory Hq is the quasifree heavy quark state and is the pole 
mass rn^^'' to this order; w(e) at e > is a perturbative probability to produce Q and 
a number of gluons in the final state which belongs to the perturbative continuum. 
It is a continuous function (plus 5(e) at e = 0). Both numerator and denominator 
in the definition of A are finite. All quantities are observables and, therefore, are 
manifestly gauge-invariant. For example, to the first order 

J pert 9 vr 9 TT 

MS/„. N / 1 , '* 16 /i 



+ . (35) 

The IR renormalon-related problems encountered in calculation of toq^'^ cancel in 
Eq. ( plf ) against the same contribution entering A(/i) in Eq. (|3^) . 

In a similar way one can define matrix elements of renormalized at the point ix 
heavy quark operators |jl^, ^ The zeroth moment of the structure function 
is proportional to {Hq\QQ\Hq) / (2Mhq) and to the short- distance renormalization 
factor. In the limit uiq oo this matrix element equals unity up to terms fi'^/niQ. 
Then it is convenient to define matrix elements of higher dimension operators in 
terms of the ratios (for which all normalization factors go away): 

. <^f<;;f';gf-'>- ^ .,,„ .„„ li^JJifM-^ . (36) 

{Hq\QQ\Hq)^ ^/-0 mQ^oo^;2 JM dey;(e, rngt;) 

_ l {HQ\QiDE)Q\HQ), _ 3 J,^ dee'wje, niQv) 

2 {Hq\QQ\Hq), - H /„-de^e,mQ^;) ' ^''^ 

and so on. All these objects are well-defined, gauge-invariant and do not depend on 
the weak current probe. 

Returning to the heavy quark mass, the literal definition given above based 
on Eg. (^41) is most suitable for /i ^ ttiq when the terms ~ {as/'n'){^'^/mQ) are 
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inessential. One can improve it including higher-order terms in fi/niQ, if necessary, 
in the following general way. 

Consider the l/mg expansion of the heavy hadron mass 



spm— av 



(38) 



We took the average hadronic mass over the heavy quark spin multiplets here. 
Besides ^1 and pf, defined above the new quantity enters the l/rriQ term. It is 

a non-local correlator of two operators Q{aiD)^Q defined in Eq. (24) of |T^. We 
apply the relation ( ^8l) to the perturbative states and get 



rri' 



pole 

Q 



A(/i) 



pert 



pert 



pert 2mQ{jJ,) 



.. . (39) 



For example, through order l/mg one has 



pert 



(40) 



To renormalize the non-local correlators entering the definition of p^, one represents 
them as a QM sum over the intermediate states and cuts it off at En > Mhq+P^- For 
example, [p^{p)]pert — (8as/97r)/i^. One subtlety is worth mentioning here: including 
higher orders in l/mq and going beyond one-loop perturbative calculations one must 
be careful calculating the coefficient functions in the expansion Eq. (|38|) . 

The general definition above is very transparent physically and most close to the 
Wilson's procedure for a theory in Minkowski space. It corresponds to considering 
the QM nonrelativistic system of slow heavy flavor hadrons and literally integrating 
out the modes with uo = E — itlq > p. Considering the SV limit and keeping only 
terms linear in v'^ is important: in this case no problem arises with choosing the 
reference energy in e or in imposing the cutoff (energy vs. invariant mass, etc.) - 
they are all equivalent in this approximation; the differences appear starting terms 
~ v^. The conceptual details will be further illuminated in subsequent publications. 



5 Determination of \Vcb\ (iKfel) from the inclusive 
widths 

We now address the central phenomenological problem - how accurately one can 
extract \Vcb\ and \Vub\ from inclusive semileptonic widths. 0{as) terms are known, 
and higher order ones are expected to be small for b decays. However radiative 
corrections to the mass get magnified due to the high power with which the latter 
enters the width. This concern has been raised in pO[. The effect of large n indeed 



shows up as the series of terms ~ (ncts/vr) . Therefore one needs to sum up such 
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terms. The resummation of the running effects in the form ^ ('^v') been 
carried out in |10[ (the term with k = 1 in this series was earher calculated in p5[|); 
with f3o/2 = 4.5 , accounting for the large-n series seems to be at least as important. 




We emphasize that the BLM-type |^6[ improvements leave out these large-n terms. 

In this paper we discussed a prescription automatically resumming the (nas/vr)^ 
terms. This eliminates the source of the perturbative uncertainty often quoted in 
the literature as jeopardizing the calculations of the widths.^ Let us demonstrate 
this assertion in a simplified setting. 

Let us consider two limiting cases, <^ m?, and ~ nib. Dropping from the 
widths inessential overall factors Vqfcp/(1927r^) ( GplVqbl"^ / {Ibn^) ) we get 

(41) 

with ai = — I ^TT^ — (or ai = —1, respectively) at = 5. As was shown before, 
these coefficients do not contain factors that scale like n, nor 02 contains n^, etc., 
and we can neglect them altogether. However, when one uses, say, the MS masses, 
one has 

^ mQ^l-^^,) (42) 
and to the same order in one arrives, instead, at a different numerical estimate 



3 TT / V 3 TT 



„ / _ {n+ l)n /4a,\^ ^ {n + l)n{n - 1) Ma, 



A is written in terms of the MS scheme masses. The expressions in Eq. (^I|) and 
Eq. (^2]) are equivalent to order a^, yet start to differ at order due to n-enhanced 
terms. The size of the difference has been taken as a numerical estimate of the impact 
of the higher order (non-BLM) corrections which then seemed to be large indeed: 
the coefficient in front of (as/vr)^ is more than 25; only such huge enhancement could 
lead to an uncertainty in the width ~ 15% from the higher order corrections. 

Since we are able to resum these terms, the width Eq. (^) has small second (and 
higher order) corrections. Alternatively, using the MS masses one is bound to have 
large higher-order perturbative coefficients, and their resummation returns one to 
using the low-scale masses. 

We, thus, conclude that the uncertainty associated with the higher-order pertur- 
bative corrections discussed in Ref. |2^ is actually absent and is an artifact of the 



^These uncertainties do not apply to our previous analyses |16, |l^ where we have employed 
quark masses evaluated at a low scale keeping in mind the natural momentum scale inherent for 
the actual inclusive widths, following the arguments under discussion. They are intrinsic, however, 
to some alternative estimates which can be found in the literature. 
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approach relying on the MS masses without a resummation of the large-n effects. 
A similar numerical uncertainty, in a somewhat softer form, appeared as a differ- 



ence between the "OS" and the "MS" all-order BLM calculations of Ref. [|TOl; it is 
likewise absent if one uses the proper low-scale masses. 

In our asymptotic treatment of the large-n limit we cannot specify the exact 
value of the normalization scale /i to be used for masses: it can equally be 0.7 GeV 
or 1.5 GeV; the exact scale would be a meaningless notion. We also saw that the 
position of the saddle point varies depending on the structure function considered. 
All these nuances are rather unimportant in practice: the impact of changing /i has 
been studied in Ref. |19|, and it was found that the values of \Vqb\ extracted from 



the widths vary by less than ±1% for any reasonable value of /i. 

The large factor of 5 also enhances the sensitivity of | V^^l to the expectation value 
of the kinetic energy operator /i^ in the currently used approach where rrib — rric is 
related to Mb, Mb*, Md, Md* and fi^. This is certainly a disadvantage; however, 
the emerging dependence of \Vcb\ on /i^ is practically identical (but differs in sign) 
to the one in the determination of |Vcf,| from the zero recoil rate B — > D*ii' (for 
details see the recent review [^). We do not see, therefore, a way to eliminate 
the uncertainty due to fi^ other than to extract its numerical value from the data. 
Having at hand its proper field-theoretic definition, one does not expect significant 
theoretical uncertainties in it. In the exclusive case, unfortunately, even this would 
not remove the sizable element of model-dependence. 

The 1/n expansion can equally well be applied when the final-state quark mass 
is practically zero as in 6 ^ u. The large corrections of order that appear when 
rrib is taken at the high scale turn out to be quite small when is evaluated at 
A/n ~ 1 GeV with A denoting the energy release. Using the analysis of Ref. fiW \ 
one can calculate the inclusive total semileptonic width r(i? — »• IvXy) quite reliably 
in terms of |\^(u6)| and mj,(lGeV) as inferred from T spectroscopy. Numerically 
one has 

. 0.00415 m'^^^J")) ' (l^) ' . (44) 
' ' V 0-0016 I \ tb ) ^ ' 

The theoretical uncertainty here is smaller than the experimental error bars which 
can be expected in the near future. 



6 Conclusions 

There are three dimensional parameters relevant for semileptonic decays of beauty 
hadrons: nib, the energy release A;,, = nib — and Aqcd- Since nib, ^bq ^ Aqcd, 

more complete BLM calculation has been carried out in the numerical results were 
quoted, however, only for the OS and MS schemes; the latter is expected to suffer from large 
rt-related corrections, whereas the OS scheme is affected by the IR renormalons when higher-order 
BLM corrections are included; the BLM resummation does not cure it completely, again due to 
the n^-enhanced terms in the spurious 1/m^ corrections. 
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the heavy quark expansion in terms of powers of Aqqb /^bq yields a meaningful 
treatment, with higher-order terms quickly fading in importance. It is natural then 
that uncertainties in the perturbative treatment limit the numerical reliability of 
the theoretical predictions. It was even suggested that the latter uncertainties are 
essentially beyond control. 

Our analysis shows that such fears ggerated. First, we have reminded the 

reader that the dependence of the total width on the fifth power of rrif, largely reflects 
the kinematics of the lepton pair phase space. We then used an expansion in n to 
resum higher order contributions of kinematical origin by identifying unequivocally 
the relevant dynamical scales at which the quark masses have to be evaluated. 
We found that the relevant scale is not set by the energy release, but is lower, 
parametrically of order Ahg/n. 

We demonstrated that the expansion in 1/n works reasonably well in simple 
examples. This is not the main virtue, in our opinion. More important is the 
fact that this expansion allows one to have a qualitative picture of different types of 
corrections based on the scaling behavior in n. In particular, it shows the emergence 
of new, essentially lower scales relevant in the semileptonic decays. 

Even without a dedicated analysis it is obvious that the typical scale of the energy 
release in the semileptonic b decays lies below rrib. Without a free parameter at hand 
it is not too convincing to defend say, the scale mb/2 as opposed to 2mb. Needless to 
say that this scale variation has a noticeable impact on the final numbers. Using n as 
an expansion parameter the ambiguity is resolved - the appropriate normalization 
scale for masses is even below mb/2 - and already this, quite weak, statement is 
extremely helpful |]19| in numerical estimates. 

Based on the large-n expansion, we arrived at a few concrete conclusions. 

• One can resum the dominant n'' terms in the perturbative expansion of the 
inclusive widths, by merely using the Euclidean low-scale quark masses (e.g. nor- 
malized at the scale ~ A/n). Therefore, the previous calculations of the width 
||T6| , are free of the large uncertainties noted in pO| and which were later claimed 



to be inherent to the inclusive widths. 

• We introduced the so called "extended" SV limit. We show that rric need 
not be close to rrify for the SV regime to emerge in the b ^ c inclusive decays. 
Large n helps ensure this regime, which gives a rationale for the relevance of the SV 
consideration in the actual inclusive b ^ c decays, both at the quark-gluon and at 
the hadronic levels. 

Recently, the complete second-order corrections to the zero-recoil formf actors 
have been computed |]2^; the non-trivial non-BLM parts proved to be small. The 



arguments based on the large-n expansion suggest then that the non-BLM pertur- 
bative corrections not computed so far for rsi(6 c) are not large either. 

The approach suggested here is applicable to other problems as well. For ex- 
ample, we anticipate that the second-order QED corrections to the muon lifetime 
which have not been calculated so far, will not show large coefficients if expressed in 
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terms on the muon mass normalized at the scale ~ vn^j?) — moreover, using 

ttem at a similar scale (in the V scheme) is also advantageous - though, clearly, it 
does not matter in practice in QED. 

Note added: When this paper was prepared for publication the full two-loop 
calculation of the perturbative correction in 6 — > c £z/ at another extreme kinematic 
point = was completed p8|. The result suggests that the non-BLM perturba- 



tive corrections to the width are indeed small if one relies on the low-scale quark 
masses, in accord with our expectations. 
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Appendix: The refined 1/n expansion 



Let us start from Eq. (0) expressed in terms of the quark masses mh and mc, with 
A = rrib — rric, and introduce variables r] and uj defined by 



A ' ^=A 



— 2mbqo + — rric 



(45) 



Thus, UJ measures the effective mass Mx in the final hadron state, uj = {Mx—rric) / A, 
(in particular, uj = represents the free quark decay) whereas rj determines the 
difference between and A^. Then we have 



r] [rj + 



7(n 
2mr 



nib 



A 





-v) 






A 


2mbJ 


Arrib 



n-2 



X 



A^ 

3ml 



1^ dr]r]^/\l-r]Y-^ (^r] + 



2m, A 

'a' 



1 - 



2mh 



vi{r]) + 
A 



Am, 



■T] 



3/2 



where functions Vi{ri) are defined through 



(46) 



vi{ri) ■ {2mcri + rj ) 



2U/2 



2mi 



2 A n 1/2 
4mh 
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nih Jo 



1 



2mh 



2771, 



■iU + 



27r 



(47) 



V2{r]) ■ {2711^7] + 7]"^ Af/^ ■ 

A n 



1 ^ I 1 ^ 7^ V + 

2mb/ 4mfe 



3/2 



THb JO 

(t] - u;)(2mc + 77A + ujA) ( 1 - ( 1 



do; (mc + ujA.) x 



2mJ ' - 



2mft 



A(?72 - u;2 
4mft 



3/2 



1^2 

2^ 



(48) 

Equation (^) represents an identity with two gratifying features. First, the 
width is expressed through the quantities fj. Being weighted integrals of the struc- 
ture functions Wi they are smoother analytically than Wi themselves. In a certain 
respect, Vi are generalizations of the moments Jj, that are relevant for calculating 
inclusive width. 

Second, the form of Eq. ( ^61) is particularly well suited for deriving the large- 
n expansion. One typically encounters integrals of the form /Q^d?7?7"(l — 7]Y . At 
k ^ 00 the integral is saturated at rjo 1 — a/{k + a) 1 — a/k, i.e. 1 — tiq <ti I. 
(We loosely refer to this saturation as a saddle-point evaluation, although it is not 
really a saddle point calculation in its standard definition.) This was actually the 
approximation used in the simple large-n expansion. However, for k = 0,2 and 
a ~ 1/2^3/2 the ansatz 770 ^ 1 is quite poor numerically. On the other hand 
Jq d?7?7"(l — 7])^ still has a reliable 'saddle point' even for large a and = 0, since 
the width of the distribution is governed by ak/ (a -|- kY and does not become large. 
The point only shifts somewhat upward as compared to the 'naive' approach. This 
is the idea lying behind the improved expansion. Its purpose is merely to determine 
the essential kinematics in the process at hand, and to use the hadronic averages 
expanded around it. Of course, the phase space integrals always can be taken 
literally for any n, if necessary. 

As it was with the simple n expansion, one finds that the explicit dependence 
on 71 is contained in the kinematical factors that can be treated separately from the 
QCD dynamics contained in the hadronic averages fi,2(^)- There is some residual 
dependence on n entering through the value of the scale rjo at which f 1,2(^7) are to 
be evaluated (i.e., the exact shape of the weight functions). Again, one has to treat 
the vector and axial current contributions separately for 6 — c and h ^ u. 

In Table |I| we compare the exact results with those obtained from the refined 
1/n expansion for n = 5 as a function of rric/mb in the simplest setting of the tree- 
level decay where the cumbersome factors in Eq. (^) are replaced by their values 
at the 'saddle' point. We have used the following notation there. The quantities 
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mc/rrib 





0.2 


0.3 


0.4 


0.6 
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-DJA,]^ n— 5 


0.223 


0.524 


0.462 


0.462 


0.476 


0.488 


0.494 


Rr)(A) 

-'-'-'^1, exact 


0.25 


0.396 


0.434 


0.459 


0.486 


0.497 


0.5 




0.229 


0.093 


0.060 


0.0374 


1.26 ■ 10-2 


2.48 ■ 10-3 





-'-'-'^l, exact 


0.25 


0.104 


0.066 


0.041 


1.36 ■ 10-2 


2.65 • 10-3 





Rr)(A,y) 


0.171 


0.248 


0.24 


0.233 


0.225 


0.219 


0.215 


■'-'■'^2, exact 


0.25 


0.25 


0.25 


0.25 


0.25 


0.25 


0.25 


Tsl, n=5 / Tsl, exact 


.80 


1.11 


1.00 


0.97 


0.94 


0.93 


0.924 



Table 1: Comparing widths in the refined n expansion with the exact results at 
n = 5] the cases rric = and rric > 0.2 rrif, are treated separately. 



Ii2^^ denote the width factors for the vector and axial contributions obtained by 



integrating wi and W2, respectively. We then have 7si = 
One defines branching ratios normalized to the exact {n 



5) tree level expression: 
(49) 



7sl, exact ("^cMfe) 

and calculates them using, on the one hand, the expansion in n evaluated at n = 5, 
and on the other hand the exact n = 5 results. Since the axial and vector part of W2 
coincide in the tree approximation, only one of them is shown in Table 0. One sees 
that the n-expansion works with a typical accuracy of about 10% for the inclusive 
width, as expected. The weight of the nonleading terms decreases as rric ^b, i-e. 
in the SV limit. It is remarkable that the expansion based on the SV kinematics 
works so well down to a rather small mass ratio of rriq/mb ~ 0.2! This illustrates 
the observation made above that it is the parameter {ruh — rric) / nrric that describes 
the proximity to the SV limit. For further analysis we note that at rric/mh = 0.3 
the 'saddle points' for n = 5 occur at r^* = 1 — \^/{mh — rric) equal to 0.16, 0.28 
and 0.4 for 7i^'*, 7i^^ and 72, respectively; for the h ^ u case r/* is 0.33 and 0.55 for 
71^2, respectively. 

Another cross-check of the numerical reliability of the refined 1/n expansion at 
= 5 is a comparison with the known perturbative expression at the one-loop level. 
In the SV regime one finds through order ag/vr 



7si,sv 



15 



a. 



0.924 - 0.945 ■ — = — ■ 0.924 1 - 1.023 • — 



TT / 15 

to be compared to the exact result to that order 



TT 



(50) 



7si .exact 



(51) 



-AM 1 - ^ 

15 V TT 

i.e., a difference of two percent only. 

From these comparisons we conclude that the refined 1/n expansion yields good 
numerical results already for the physical value n = 5. 
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